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I , guirmw.-ry 


The purpose of this contract was to expand the epicyclic gear dynamics 
program to add the option of evaluating the tooth pair dynamics for two 
epicyclic gear stages with peripheral components. The option was 
developed for either stage to be a basic planetary, star, single exter- 
nal-external mesh, or single external-internal mesh. The two stage 
system includes an input mass and shaft, an output mass and shaft, and 
a connecting shaft, where the shafts are each modeled with torsional 
springs and dampers. The solution procedure was nearly the same as the 
procedure previously used for determining tooth pair displacements ani 
stresses in single stages. The primary differences were that the indi- 
vidual gear component displacements were calculated rather than the net 
sun-planet or ring-planet displacements. This was necessary in order 
to determine the relative displacements between the shafts and the 
input and output gears. This generally increases the number of degrees 
of freedom to be solved per stage compared to the single stage solu- 
tion. 

The option to evaluate two stages makes the user's job more difficult. 
Two stages of basic gear system information must be input as well as 
the additional component information. In addition, the boundary condi- 
tions and associated iteration procedure become more complex. Thi« is 
due to both the i n creased number of components and to the ti me for a 
complete mesh generally being different for each stage. 

A brief investigation into methods of reducing the program's computa- 
tion time was done. The efforts focused on reducing the munh^p of 
iterations required for boundary condition convergence. It was recom- 
mended that the beginning and ending values be weighted differently to 
u t ilize previous iterations more effectively . As over 90% of the com- 
putation time is in the numerical integration routines, reducing the 
number of solution integrations would yield a direct, Hrvm-r reduction 
in time. 

Execution of the initial test case indicated an instability in the 
solution. The tooth pair load pattern is reasonable; however the mag- 
nitudes of the tooth pair loads grow to excessive, unrealistic values 
as a function of time. This could be due to the irvi t.-tAi conditions, a 
code error, or some type of numerical instability. A procedure has 
been recommended for eliminating possibilities and determining where 
the problem lies. 
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Over the past several years, NASA has developed a gear dynamic analysis 
arrf computer code for standard and high oontact ratio gears. Hie ana 
lvsis was expanded to i nclude internal involute tooth forms in addition 
to external involute tooth forms, then expanded to i n cl u de several pla- 
nets in epioydio gear systems. The development continued with the 
addition of 1 rwi and double h el ica l gears, a floating sun gear 

option, a natural frequency evaluation, a refined helical gear ocra- 
pn*-iw» routine, and a flexible oarrier evaluation. 

The ijoi tiA-i program was developed for a single spur gear mesh and to 
operate over a wide range of oontact ratios (up to 4.0) for analysis of 
high oontact ratio as well as low oontact ratio gearing. This single 
mesh program was an extension, by Cornell and Westervelt (references 1 
and 2), of the basic concept developed by Richardson in 1958 (reference 
3) The tooth pair compliance and stress sensitivity formulation of the 
sir^le spur gear mesh program was used in the epicyclic gear dynamics 
program, applied to each mesh. 

'TMg oontract, NAS3-25281, developed an option for two spur gear stages 
with input and output shafts with attached masses and a shaft oonneot- 
incf the two stages. This was a practical extension to the program as 
more than one g ear ing stage are often used for speed reduction, space, 
weight, and/or auxiliar y units. Thus, this extension allows for model- 
ing of peripheral components in order to assess the impact on the gear 
tooth stressing. Although the basic oode was completed, the oause of 
the unexpected instabili ty in the results oould not be deter m i n ed with 
the cost constraints of the oontract. NASA, therefore, decided to 
aooept deliv ery of the code and documentation rather than extend the 
contract at this ti m e. 



III. Task I-Multlple Gearing Stages and System Components 


The first two tasks of contract NAS3-25821 were to modify the existing 
Multi -Mesh Gear Dynamic Analysis Code (GEAREKNMULT in NASA notation, 
F178 in Hamilton Standard notation) to include a second stage and per- 
ipheral components. These tasks included the technical development and 
most of the coding efforts. Another smaller task was to investigate 
methods of reducing computation time. 


A. Analytical Model 


The dynamic response equations for a single gear stage were expanded to 
I n cl u de additional degrees of freedom for an input shaft and attached 
input mass, and an output Shaft with an output mass. Thi.q was accom- 
plished by adding two new equations and modifying other equations. Fig- 
ure 1 shows a diagram of the two stage system with no details for the 
individual gear stages. Figure 2 shows a general schematic for either 
epicyclic gear stage. The specifics of the gear tooth dynamics and 
program capabilities for each stage are discussed in References 1-7. 

For the numerical model, the input and output shafts each required an 
ad d itio n al second order differential equation to account for the iner- 
tias on the input and output shaf ts - 

The equation for the input shaft and attached mass is: 


Cl) 


where 6 inl = , Y lnpnt gfta.T fo r stage 1 

base radius of the input gear to stage 1 


The equation for the output shaft and attached mass is: 


^out ®outsh + ^out^outsh ®out2^ 

+ *out ^ 9 outsh ~ 0 out2^ “ T out 


( 2 ) 


The addition of these equations changes the way the torque is fed into 
and out of the individual stages. Previously, there was a constant 
input torque term in the sun gear equation and a constant output torque 
term in either the ring or carrier equation. Thp* torque is now fed 
into the input gear equations by way of shaft rotational displacements 
and velocities, which vary with time. The rotational displacements are 
transformed into equivalent displacements along the line of action for 
each mesh of the input gear. The output gear equations also have the 


3 



equivalent shaft displacement and velocity terms in place of the output 
torque terms. 

The new two stage program has the capability of solving for the dynamic 
tooth loads in planetary systems, star systems , single mesh external- 
external involute tooth forms, and/or single mesh external-internal 
involute systems. Tabl e 1 summarizes the a v a il a ble spur gear systems 
for two stage dynamic solutions. The two stage program has been 
expanded to allow the user to specify the input and output gears. The 
single stage solution assumed the systems were speed reduction systems, 
but the two stage solution also allows for speed increasing systems. 
Th-t.g option will allow the user to model a back-to-back test rig, which 
is generally one reduction system and a second simil a r system which 
increases the speed. Tor example, a single stage planetary system nor- 
mally has the sun gear as the input and the planet carrier for the out- 
put. For the two stage solution, the planetary system oan have the 
input component as either the sun gear or the carrier with the output 
being either the carrier or the sun gear respectively. 

Th«» equations for each separate gear component were the sa m e as pre- 
viously developed. The basic equation development was documented in 
References 4-7, and the equations are summarized in App endix A for 
convenience. The input and output gears for each stage will depend on 
the specified system type, see Table 1. Thus, the additional terms 
neerigrf in the gear equations, due to the shafting, are written in terms 
of general input and output gear equations, and are as follows. 

input gear equation for stage 1 
+ ^in ^inl - ^insh^ + ^in ^inl ^instP 


output gear equation for stage 1 
+ C 12 ^outl ~ ^±n2^ + k 12 ^ Y outl ~ Y ln2^ 


input gear equation for stage 2 
+ C 12 ^in2 - ^outl'* + k 12 ^in2 - ^outl^ 


output gear equation for stage 2 
+ ^out ^out2 - ^outsh^ + k out ^out2 ^outsh^ ^ 


where. input gear to stage 1^ ®insh 

^outsh ” output gear from stage 2^ ^outsh 


(3) 


(4) 


- 0 


(5) 
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Transformations from shaft rotations to displacements along the appro- 
priate line of action are necessary to determine the tooth pair dis- 
placements along the line of action. This allows the tooth pair dis- 
placements and loads to be determined by a procedure si m ila r to the 
single stage procedure. 

The numerical solution of the dynamic equations utilizes IMSL (Interna- 
tional Mathematic and Scientific Library) numerical integration rou- 
tines. These routines solve a system of first order differential equa- 
tions. Appendix B shows the reduction of the second order equations to 
first order equations. 


The tooth pair displac ements along the line of action are calculated 
from the gear component displacements via: 


?sp = Y s ~ *c - V w 

^rp = ^pi - Y c 7r 

For the single stage solution, the generalized coordinates to be solved 
by the numerical solution were the tooth pair displacements along the 
Htv> of action (y and y ), which minimized the number of equations 

to be numerical l.y solved. For the two stage system of equations, the 
individual gear component equations must be solved. This was neoessary 
in order to obtain the input and output gear displacements (y , y , y , 


Yp^) for both stages. Therefore the relative displacement terms bet- 
ween the shafts and the input or output gears oould be calculated. 
This means there are more degrees of freedom to be solved for each 
stage in the two stage solution . There are (3 + number of planets) 
degrees of freedom for each stage plus two degrees of freedom for the 
input and output shafts. Thus, there are a total of (8 + + N^) 


second order equations and 2(8 + N. + N~) first order equations to be 
solved simultaneously. 


The input and output shaft damping terms are calculated from a damping 
ratio, stiffness, and inertias. The damping coefficients are then cal- 
culated using the following formulations. 


c i* - 2 5 [>W < l/J m + ] 1/2 

°out - 2 ? IW ( 1/J out + "W I 1 ' 2 

c i2 - 2 ? t k i2 / < + 1/J na > ) 1/2 


(9) 

( 10 ) 

( 11 ) 
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B. Solution Prooedure 


The solution prooedure used to determine the two stage dynamic gear 
tooth loads and stresses was very similar to the solution for the 
single stage dynamic gear tooth loads and stresses. To minimize the 
program modifications, as many of the existing methods and procedures 
as possible were used. Most sections of the program needed some modi- 
fication in order to process the additional stage. Many of the 
slightly modified routines were used twioe, onoe for each stage, by 
sending the necessary information for the relevant stage. 

The single stage solution prooedure will be reviewed first, beoause the 
two stage solution prooedure is essentially an extension of the single 
stage prooedure. The general program prooedure is illustrated via a 
flowchart in Figure 3. In brief, the program calculates tooth pair 
loads and solves the differential equations for each of 100 time 
steps — a pieoewise linear solution. The boundary conditions are deter- 
mined via an iterative prooedure, which will be explained in more 
detail later. 

The boundary condition iteration prooedure is intended to lead to a 
steady state condition. For the single mesh this is accomplished by 
evaluating one full mesh cycle, where a mesh cycle means the time 
required for a tooth pair to move the length of the of action. 

The mesh cycle time for spur gears is divided into 100 time steps for 
the numerical time stepping solution, where each increment of time is 
treated as linear with respect to the tooth pair stiffnesses. The 
boundary conditions, for sun-planet and ring-planet tooth pair dis- 
placements and velocities, are compared at the beginning and end of the 
mesh time cycle. If the displacements and velocities are within a spe- 
cified toleranoe, the system has reached a steady state. If the condi- 
tions are not within the toleranoe the initial and final values of the 
displace m e nt s and velocities are averaged to obtain a new set of ini- 
tial conditions, and the prooedure Is repeated with new ini tial condi- 
tions. This prooedure is repeated until the boundary oaniitions are 
all within the specified toleranoe. The final iterated 100 time step 
solution for tooth pair loads and displacements is then post-prooessed 
for tooth stressing. 

The boundary conditions for two stages become quite complex compared to 
one stage. The main complexity is due to the two stages having differ- 
ent mesh time cycles. This means that steady state is no longer 
achieved for one mesh cycle, rather it must he a mutually common timp 
for both stages, analogous to a common denominator. Figure 4 schemati- 
cally shows two examples. The first one shows a case where the two 
stages have time cycles that oan readily achieve a steady state, i.e. 
one stage has a time cycle exactly 2 times the other stage's time 
cycle. The other example Shows a oase where 47 cycles would have to be 
evaluated for steady state, which would consume large quantities of 
computer time (and a code change to increase the size of an array). 
However, it is thought that this oase could be approximated by evaluat- 
ing five cycles. It follows that there are systems in which there is 
no true steady state. 
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Qnoe the boundary conditions have converged, the tooth pair steady 
state displacements are post-prooessed for stressing through the 100 
time steps. If any tooth pair errors have been requested for evalua- 
tion, the program continues to evaluate 10 time cycles to simulate the 
user input tooth pair errors going through a mesh for either one or two 
stages. The method used to prooess the tooth pair errors starts with 
the converged boundary conditions from the no tooth errors solution as 
the initial conditions for the numerical solution with errors. The pro- 
gram thus simulates a transient response (for a duration of ten mesh 
cycles) as the various tooth pair errors came into mesh. The dynamic 
effects of the tooth pair errors have generally been observed to dampen 
out after evaluation of 3-5 mesh cycles. 

For the two stage expansion, it was desirable to work with the current 
procedure to maintain commonality. This implied that the 100 time 
steps, corresponding to a mesh cycle time, needed to be maintained, as 
this fundamental approach is used throughout the program. The first 
question was which stage's time cycle should be used to correspond to 
the 100 time steps embedded in the program. A smaller time step for 
the numerical solution would tend to make for a more stable numerical 
solution, while the larger time step would reduce the computation time. 
The smaller time step was chosen for stability reasons. 

The other question was how to handle the two different mesh tingy? in 
order to achieve a "steady state" in a practical amount of computation 
time. It was decided that the program could be set up to evaluate the 
boundary conditions after a user specified number of time cycles was 
evaluated. However, the fact that there could be nacaag with no tr ue 
steady state made it desirable to try and approximate a st eady state 
and then evaluate a transient response using the approximated bo undar y 
conditions as initial conditions. It is possible that the instability 
could be related to these approximate steady state boundary conditions. 

In order to verify the approximate steady state boundary condition con- 
cept, the original single mesh program was used to determine if a 
change in the mesh time cycle would significantly affect the converged 
boundary condition values. Several cases were run with slightly dif- 
ferent numbers of teeth (which determines the mesh time) . The iterated 
results showed minimal variation in the converged values for displace- 
ment and velocity. Thus, this implied that it would be po ssible to 
evaluate 'approximate two stage boundary conditions ' at a time when the 
two stages were close to a steady state. The program will use the ini- 
tial conditions obtained from the approximate steady state solution and 
continue to evaluate 10 additional cycles for a transient response. 

The transient response utilizes the same procedure used to evaluate 
tooth pair errors in the single stage code. The solution procedure, 
with 100 time steps per cycle, is essentially the same except there is 
no boundary condition averaging. It should be noted that the two stage 
solution compares the tooth pair displacements and velocities for con- 
vergence, not the individual gear components. The ini t.ia.1 and final 
conditions for the individual components are each averaged and tha t 
average is used for the initial conditions for the next iteration. 


7 



The subrou tines written for the two stage solution generally follow 
similar notation to that used in previous codes, both in subroutine 
names arri. variable names. The code was written assuming a maximum of 10 
planets per stage, spur gears only, and no floating sun or flexible 
carrier option. 

Several new routines and some slightly modified routi n e s were written 
in order to accommodate the additional information required for two 
different stages. The number of planets per stage is limited to 
10, which should be more than sufficient for any real system. This is 
a reduction in the number of planets allowed in the single stage 
solution but simplif ies the oode changes by utilizing the existing 
arrays and the corresponding dime n sions. The arrays now contain stage 
1 information in the first 10 elements and stage 2 information in the 
second 10 elements, for the same total number of array elements as the 
previous solution. 

The procedures and calculations of the new program are essentially the 
saiae as the single stage program. The primary changes come through 
nearly every section or subroutine being executed for stage 1 and then 
for stage 2. Same of the basic parameters, such as mesh time cycle, 
required additional logic in order to implement the option of 
"reversed" systems (where the systems are speed increasers rather than 
reducers). Most of the minor modifications were made to accommodate 
the additional parameters of the second stage, and did not change any 
basic concepts. Beoause of the option of multiple cycles being evalu- 
ated for bo undar y ooidition convergence, it was also neoessary to add a 
dummy time parameter to obtain proper tooth pair contact . Without the 
dummy time parameter, the stage with the shorter mesh time cycle lost 
an tooth pair contact during the second 100 time steps . Hms, to 
obtain proper contact evaluation, a dummy ti m e was reset to zero at the 
beginning of each wrath time cycle for the stage with the smaller mesh 
time cycle. 

The in put routines were modified to meet NASA programming sta nd a r d s . 

The additional stage of information neoessary for the two stage solu- 
tion required modification of the code to set up the Input variables. 
This portion of the oode was previously in the preprocessing routine; 
however, the preprocessing routine (READ2) is a very large routine and 
does not meet the 150 executable line requirement. Thus, the section 
of the code for setting input variables to meaningful notation was 
modified to be separate subroutines, accepting input for either one or 
two stages of information, and the remainder of KEAD2 remained 
unchanged except for additional arguments. Three subroutines were 
neoessary to convert a section of the old routine into acceptable 
length rou tines ani a fourth was neoessary to control which stage was 
being preprocessed. 
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C. Methods of Reducing Computation Time 

Part of Task I was to recommend a method for possi b le future develop- 
ment work in "alternate methods for expanding the model to more than 
two stages (such as a superposition dynamic model) . " Hie stated reason 
for s uch a recommendation was the potential for extensive computational 
times. 

r php approach to this task was to emphasize the reduction in computer 
time required by the program. In analyzing the present CPU usage of 
the program, it was determined that over 90 percent of the computa- 
tional time was in the numerical integration routine. Hie basic alter- 
natives to reducing program solution time were as follows. 

a. Reduce the number of time steps used through the 
mesh for a direct, linear reduction in time. 

b. A more advanced solver routine for an unknown time 
reduction. 

c. Reduce the number of iterations required for 
boundary condition convergence for a direct li n ear 
reduction in time. 

Item a. was thought to be practical for single mesh oases only, since 
only in these cases could the timing of events through the mesh time 
cycle, be accurately estimated. Item b. was not thoroughly explored, 
but was felt to hold some promise. Item c. was the main focus of this 
study, ard could be applicable to any level of the oode, i.e. , single 
mesh, multiple mesh, and multiple stages. 

In order to investigate better methods of convergence, the flexible 
carrier test case given in Reference 4 was used. This case was ini- 
tially solved by calculating input conditions as specified in Appendix 
A of Reference 4, and many boundary condition Iterations, which vise 
significant CPU time. 

Figure 5 is a plot of the ending values of the sun gear dis p lacement 
a rd velocity versus the number of iterations through a gear mesh 
period. It Shows that the method of selecting new values is creating a 
cyclic behavior that will take many iterations to converge. The mean 
value of the plotted functions match the converged value of the solu- 
tion. The behavior of the planet and ring gear results was similar. 

In order to evaluate the convergence difficulties, the code was tempo- 
rarily changed to weigh the beginning and ending values differently in 
order to obtain a higher weighting to previous iterations. The revised 
code. Shown in Figure 6, has not been incorporated in the new program, 
but was vised for this investigation only. The convergence could be 
tailored for a particular case by changing coefficients of the dis- 
placements from the previous iteration and the current iteration (e.g. 
SXSO and XSO respectively), and similarly for the velocities (SCSI and 
XS1). The sum of the coefficients must be equal to unity. The results 
of using 0.7 and 0.3 for the displacement, and 0.5 and 0.5 for the 
velocity are shown in Figures 7 and 8. 
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In order to reduce the computation time for a dynamic solution, it is 
recommended that this method he added to the code, with the weighting 
coefficients as input parameters rather than fixed values in the oode. 
It is further reoammended that oode be written to automate the plots as 
shown in Figures 7 and 8. The resulting oode should substantially 
enhance the user's understanding and the oonvergenoe of the boundary 
conditions. 


D. Discussion of Instability and Test Cases 

Task Ic was to evaluate the two stage model and determine the interac- 
tion between two stages and the effects of input and output shafting 
and attached masses. This was not fully accomplished due to the ini- 
tial test case showing an instability in the results. However, same 
general comments will be made regarding the apparent instability and 
other observations from a variety of test oases. 

The first test oase that was used approximated the CH-58 dropped tooth 
planetary design being tested at NASA Lewis. A dropped tooth design 
“drops" a tooth from each of the planets so that the number of teeth on 
a planet divided by the number of planets is not an integer number. 

The purpose of the design is to reduce dynamic tooth loading. The 
actual pl ane ta r y system could not be modeled with the program, because 
the dropped tooth design leads to different pressure angles for the 
sun-planet and ring-planet meshes. The program equal pressure 

angles for the sun-planet and ring-planet mesh, although an equivalent 
buttress tooth form might be a feasible modeling method. Figure 9a 
shows the input test oase for a single stage and Figure 9b shows the 
corresponding output summary. This oase had ftflawwi the pressure angle 
of the sun-pla n et mesh for both the sun-planet and ring-planet meshra 
and used the program's geometry preprocessor. 

The next step was to evaluate two stages, sirmiiAtin g the bacik-to-bacik 
test rig at NASA. An example input data set for two stages is shown in 
Figure 10. The first stage was modeled as a planetary, with the sun 
gear as the input gear, and the second stage was a “reversed planetary" 
with the carrier as the input. This oase was used for initial verifi- 
cation of the code interaction and debugging purposes. The geometry for 
the two stages was the same, with the primary difference between the 
stages being the input torque and rptn. The different torques led to 
slightly different tooth pair stiffnesses, as the torque level influen- 
ces the Hertzian oomponent of the stiffness function, see reference 2. 
The difference in rpm oould affeot the dynamio load levels, particu- 
larly if near a resonant speed. 

Table 3 summarizes the test oases run with the two stage program and 
the resulting component displacements. The initial test oase, oa s e 1, 
was designed to uncouple stage 1 from stage 2. Thus, the connecting 
shaft stiffness was very low. The input and output shaft stiff ness es 
were also low, because some preliminary test oase results indicated 
large displacements for very stiff input and output shafts. The solu- 
tion seemed to be very sensitive to the values of the shaft stiff- 
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nesses, as indicated by comparison of case 4 and 5 (stage 1). For a 
decrease in input shaft stiffness by a factor of ten, the magnitudes of 
the planetary component displacements along the line of action 
decreased by orders of magnitude. 

Some of the dynamic load results for the initial test oase have been 
plotted in Figure 11. These results illustrate a typical tooth pair 
meshing pattern of engagement and disengagement of the teeth. The 
meshing pattern also repeats for subsequent cycles, as would be 
expected. What is unusual about these test results, is the subs tantial 
incr ease in the magnitude of the loads as time progresses. This 
I n creasi n g magnitude is also illustrated in Figure 12, which shows the 
maximum stress in a cycle vs. the transient cycle number. It is evi- 
dent from this figure that there is some type of instability in the 
solution. It is interesting that the Hertz stress increases in a 
nearly linear fashion, while the maximum bending stress is a h-ighg-r 
order function. This is logical in that a nonlinear increase in loads 
would correspond to a nonlinear increase in bending stress. It is 
recommended that the system equations be evaluated for stability. 

The cases summarized in Table 3 show the effect of varying the input 
and output masses and shaft stiffnesses. Only one boundary condi tion 
iteration was performed, which was sufficient to observe the change in 
the component displacements due to the modified pararoeter(s) . The 
affect of varying the masses on the input or output shaft was inversely 
proportional to the input or output shaft displacements respectively. 
For example, the input mass for oase 2 is 10 ttmap the input magg for 
case 1 and the input Shaft rotational displacement decreased by the 
same factor of 10. The other stage 1 component displacements, sun, 
carrier, and planets, also Showed decreases of nearly the same factor. 
This indicates that the Shaft displacements are dominating the individ- 
ual input and output gear equations, and the tooth pair dynamics are 
therefore overshadowed. 

Case 3 is the same as Case 1 except the input and output masses have 
been reversed. The input and output Shaft displacements were corre- 
spondingly interchanged. As would be expected, the gear displacements 
for the two stages were interchanged also, but with some variation due 
to the change in torque and rpm of the two stages. 

In Case 4, a very stiff Shaft connecting the two stages was investi- 
gated. This case verified that the displacement of the output g ear 
from Stage 1 (planet carrier) would equal the displacement of the input 
gear to Stage 2 (planet carrier) for a rigid connection. 

Case 5 was a check case to verify that the input tihcm Shaft cnn~M 
be isolated from Stage 1 by using zero torsional Shaft stiffness. By 
using a low connecting shaft stiffness, this case also verified that 
the two stages could be isolated from each other. By having zero tor- 
sional input Shaft stiffness, there is no torque transmitted to the 
input gear of Stage 1 and therefore approximately zero displacements in 
Stage 1. The second stage displacements were nearly the same as Case 
4, which verified the stages were isolated from each other. 
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Case 6 was nearly the same as Case 4, bat the connecting shaft stiff- 
ness was a lower and more realistic value. All of the component dis- 
placements were slightly larger for the softer shaft. There was also a 
gmai i difference in the magnitudes of the output gear from Stage 1 and 
the input gear to Stage 1. This is logical, because as the connecting 
shaft stiffness decreases, the relative displacement of the two ends 
will increase. 

Case 7 increased the torsi onal stiffness of the input shaft by several 
orders of magnitude. The resulting displacements shown In Table 4 
indicate the dis placement s also increased by orders of magnitude for 
the first stage components. The second stage component displacements 
also increased substantially, but not as much as the first stage. 

There are several possibili ties for the souroe of the Instability. 
Figure 13 illustrates a procedure for eliminating po ss i bili ties and 
tracking down the problem. The recommended procedure would start with 
a very simple "two stage" system, as shown schematically in Figure 14. 
TMg simple system, effectively a gear train, with constant tooth pair 
stiffnesses could be evaluated using a small program, separate from the 
gear program . Hie system should have a 1 • 1 gear ratio for simplicity 
to ens ure a steady state oan be achieved for boundary oondltlon 
oonvergenoe. If the simple system shows an unstable behavior, it will 
be known that there is a problem with the basic problem formulation. 
That is the least likely souroe of the problem, but it is the logioal 
starting point. The model complexity oan then be increased by adding 
in the variable tooth pair stiffness. If the independent solution 
yields r easonabl e results, the same system oan be run in the full mul- 
tiple mesh gear program. If the results are inconsistent, it would 
indicate an error in the two stage multiple mesh oode. 

If the results are consistent for the 1:1 gear ratio oases, the pro- 
Han is most likely related to the initial conditions and/or the boun- 
dary condi tion iteration scheme. This would indloate a better method 
for es timating the initial conditions and for determi n i n g oonvergenoe 
fiVimii h be devised. It should be noted that this is not minor task. In 
general, each stage plus an input or output shaft will have 2 times (4 
+ number of planets) for boundary conditions to be determined. Not 
only is it dif ficult to obtain initial estimates for the tooth dis- 
placements and velocities for each component of the epicyclic stage, 
but it is even more difficult to determine a general method of combin- 
ing the two stages. 

Another test case was tried to determine if the instability was system 
dependant. The second oase chosen was the Stoeckicht 2 K-H planetary 
from Reference 5. The two stage test oase was set up as two planetary 
reduction stages with identical geometry and basic gear parameters. 

•jhe only difference between these two stages would be the torcjue and 
the rjm. , t , m « case exhibited similar instabilities, with tooth loads 
increasing substantially with time. 
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IV. QQDQlUdiD g Remarks 


An option to evaluate two gear stages, with a connecting shaft and 
input ani output shafts with attached masses, has been added to the 
epicyclic gear tooth dynamio analysis program. There were extensive 
code modifications to accommodate two stages of: geometry preprocess- 

ing, tooth pair load calculations, and dynamic equations. The new code 
u tiliz es the same general methods and procedures as were previously 
used for the single stage solution and includes additional degrees of 
freedom for the attached masses. The two stage option also allows for 
speed increasing systems as well as speed reduction systems. 

It is recommended that code changes be incorporated for a more robust 
boundary condition iteration scheme. The proposed code changes can 
signifioantly reduoe the number of iterations required for a converged 
steady state solution, and therefore reduoe the computation time for a 
dynamic solution. 

The program's plotting capabilities should be modified to use more uni- 
versal plotting routines. The plotting capabilities should also be 
expanded to prooess the large amount of output from a two stage solu- 
tion. In addition, plotting capabilities for the boundary conditions 
should be added to enhanoe the user's understanding of the boundary 
condition convergence prooess. 

The initial test case indicated an instability in the dynamio solution, 
where the gear tooth displacements continually increase with time, bey- 
ond any reasonable values. Several variations on shaft stiffnesses, 
attached masses and shaft damping were evaluated to investigate the 
problem. The tooth pair loads with respect to time indicated the tooth 
pair contact patterns were reasonable, in that the basic pattern was 
repeated with time, although the magnitudes continued to increase. 
Further evaluation of the problem could not be made because of fund 
limitations . 
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Appendix A. Epigycllc G ear System Equations for Either Stage 


The following second order equations of motion describe the 
individual gear components in an epicyclic system. They are 
written in terms of the displacements along the lines of action. 
The development of the basic equations is documented primarily in 
Reference 1 and 6. 


Sun gear equation: 


-a + JL + ii v 


(A. 1) 


Planet gear equation: 


m 


p i\ ds P j s P i + d rp/rp ± L S p i + L rp ± 0 


Ring gear equation: 


.. N . N 

my - Id y -il 

P r i-1 rp i Yr Pi i=i r P ± 


/R b 


out 


(A. 2) 


(A. 3) 


Carrier equation: 


.. N . N , N N 


T out /R b 

o 


NOTE. For the two stage solution, the constant torque terms irrii- 

g ° t0 Z ^ ro - ^ is input to the appropriate 

gear equations via terms for relative displacement with respect 

tions° 0th Pa±r l0ad termS 8X16 determined from the following equa- 


m 


1 - e ~ - X 2 - S 2 _ )u ® J (A. 5) 


SPi jti 'sPi sp 


JI 


spji spji spji 


m 


L r Pi ‘ - % - *%> 6% % "Vp 1 (A - 6 > 


Ji Ji ji *ji 
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Appendix B: First Order Differen tial Equations 


The numerical integration routines solve for a system of first 
order equations. The transformation of the second order differ- 
ential equations to first order equations is shown below. The 
subscripted numbers on the derivatives correspond to elements in 
the array (XI) in the program. 

Input shaft equation reduction: 


*3 " W 

*4 - *3 * S insh 
Then , , 

*5 - t V - <W*4 - W " k l* (X 4 - 

*4 ‘ *3 


(B.l) 

(B.2) 

(B.3) 

(B.4) 


Output shaft equation reduction: 

Let , _ .. 

x 5 “ 9 outsh 
Xg = x g = e 0UtstL 

Then t 

*5 - 1 T out " C OUt (X 6 

^6 “ X 5 

The remaining equations will vary depending on the specific sys- 
tem type, see table of system types in the main report text. For 
simplicity, the following first order equation derivations will 
aggnmp both stage 1 and stage 2 axe planetary systems. Thus, the 
generalized input aid output gear equations in equations 3-6 
are now specified, where the input gear is the sun gear and the 
output gear is the planet carrier for both stages. 

Sun gear equation for stage 1: 

Let 

x l " Ys " ®inl *be 


(B.9) 
(B. 10) 


(B. 5) 
(B.6) 

- 6 out2 ) " k out Cx 6 “ 0 out2^ /J out (B - 7) 

(B.8) 
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Then 


N1 N1 



+ C lA - W 



- k i!> Cx 2 - "reV I ' m s 

(B. 11) 


*2 ‘ *1 

(B. 12) 

Carrier equation for stage 1: 


Let 

« • • 

^ “Ye 

(B. 13) 


*8 = *7 = 

(B. 14) 

Then 

N1 N1 N1 N1 

x_, = [ £ d w y w + I d_ y + £ L_ + £ L_ 

^ i-1 ^i rp i i-1 ^i ^i i=l ^i i=l ^i 



“ C 12 ( ' X 8 “ X 12^ ~ k 12^ X 8 “ X 12^ 1 m c 

(B. 15) 


X 8 =X 7 

(B. 16) 

Ring 

gear equation for stage 1: 



* •* _ __ 

*9 - y r - o.o 

(B. 17) 


x 10 ' ^ ' % " 00 

(B. 18) 

Note: ring gear is fixed for a planetary system. 


Planet gear equations for stage 1: 


Let 

*1+16 ■ ?pi for i - 1 to N1 

(B. 19) 


X i+Nl+16 ‘ x l+16 fOT 1 - 1 t0 1,1 

(B.20) 

Then 

vie - [ vv VV V V 7 \ 

(B.21) 


• 

X i+Nl+16 = x i+16 

(B.22) 
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Sun gear equation for stage 2: 
Let 



X 11 

II 

to 

to 

(B.23) 


X 12 

“ X 11 - Y s2 

(B.24) 

Then 

X 11 

N2 . N2 . 

- - Y sp 2i + ° 12(Xl2 " ^ 




+ k 12 (x 12 - Xg) ] / m g2 

(B. 25) 


X 12 

“ x ll 

(B.26) 

Carrier equation for stage 2: 


Let 

X 13 

» • •• 

Y c2 - ®out2 ^bc2 

(B.27) 


X 14 

= X 13 = Y c2 

(B.28) 


Then 


N2 N2 

X 13 = ^ 1 d rp2 i Y rp2 i + i ^ 1 d sp2 i Y sp2 i + 
" C out ( ‘ R bc2 X 6 ■ X 14' ) “ k out ( ' R bo2 X 6 " 


N2 

E L i 0 + £ L * 

sp2 ± ± _ 1 rp2 ± 


N2 
£ : 
i=l 


x 14 ) ] / m 


C2 


(B.29) 


14 


= x 


13 


(B.30) 


Ring gear equation for stage 2: 


x 


15 


x 


16 


Note: 


= 0.0 

■ X 1S ■ *r2 - 0 0 

ring gear is fixed for a planetary system. 


(B.31) 
(B. 32) 
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Planet gear equations for stage 2: 
Let 


X i+2N1+16 X pi 


for i = 1 to N2 


X i+2N1+N2+16 X i+2N1+16 


X i+2N1+16 I d sp2 i Y sp2 i d rp2 i Y rp2 i + L sp2 i 


- L_ 0 ] / m 0 
rp2 ± J p2 ± 


X i+2N1+N2+16 


X i+2N1+16 


(B.33) 

(B.34) 


(B.35) 

(B.36) 
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Appendix C. 


= damping coefficient for the input shaft (in-lbs-s/rad) 

C Qut = damping coefficient for the output shaft (in-lbs-s/rad) 

C = damping coefficient for the output shaft (in-lbs-s/rad) 

d d = tooth pair damping (lb-s/in) 

rp: sp. 

e , e = tooth spacing error (in) 

r? ij **ij 

J = rotational inertia attached to input shaft (in-lbs-s 1 /rad) 
in 

J = rotational inertia attached to output shaft (in-lbs-s Vrad) 
out 

k = torsional stiffness of the input shaft (in-lb /rad) 
in 

k = torsional stiffness of the output shaft (in-lb/ rad) 

‘out 

k = torsional stiffness of the shaft connecting the two stages 
X 2 

t l tooth pair loads for planet mesh i (lbs) 

r Pi' ^i 

m = rotational (equivalent) mass (lb-s 1 /in) 

N = number of planets 

R = base radius (in) 

y = displacement along the line of action (in) 




cam modification, ring-planet mesh 
cam modification, sun-planet mesh 
ring-planet tooth pair spring rates 



sun-planet tooth pair spring rates 


0 = rotational displacement (rad) 
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jf = damping ratio 


= input torque (in-lb) 

T out = ou ’t'P u " t torque (in-lb) 

0 rt) = identity function for ring-planet tooth pair contact 
<£> = identity function for sun-planet tooth pair contact 

* ji 

= 0 or 1 depending on whether the tooth contact is on 

the profile modification cam or not 

Subscripts 

insh = input shaft 

outsh = out put shaft 

ini = input gear to stage 1 

in2 = input ge an to stage 2 

outl = output shaft from stage 1 

out2 = output shaft from stage 2 

in = Regarding input 

out = Regarding outp ut shaft 

s = sun gear 

pi = planet gear i 

r = ring gear 

c = carrier 

sp = sun-planet mesh 

rp = ring-planet mesh 

12 = shaft connecting stage 1 to stage 2 

1 = stage 1 

2 = stage 2 
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Table l: Spar r ’ ftar system s Ava il a b l e for Two Stage Dynam i c 

Solution 


System 

nil toV^T* 

Spur Gear 
Svstem Type 

Tnput Gear 

Output Gear 

IlUIUUCl 

1 

planetary system 

sun gear 

planet carri< 

-1 

planetary system 

planet carrier 

sun gear 

2 

star system 

sun gear 

ring gear 

-2 

star system 

ring gear 

sun gear 

4 

single external- 
external* mesh 

sun gear 

planet gear 

-4 

single external- 
external* mesh 

planet gear 

sun gear 

5 

single external- 
internal** mesh 

planet gear 

ring gear 

-5 

single external- 
internal** mesh 

ring gear 

planet gear 


* External-external mesh means a pinion and a gear, both with 
external involute tooth forms. 

** External-internal mesh means a pinion and a gear, with the 
pinion having an external involute tooth form and the gear 
having an internal tooth form. 
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Table 2 : OH-58 Planetary Characteristics 


System Type 
Speed Range 
Torque Range 
Number of Planets 


Planetary 

Up to 1620 rpm (sun gear) 

Up to 12,450 in-lb (sun gear) 
3 


Stage 1 input : 
Stage 2 input : 
Stage 1 output: 
Stage 2 output: 


Sun gear 
Planet Carrier 
Planet Carrier 
Sun gear 


Peripheral 

Input 

Output 

Connecting 

Characteristic 

Shaft 

Shaft 

Shaft 

Inertia (in-lb-s 2 ) 

2.5 

1.0 


Stiffness (ln-lb/rad) 

l.OxlO 5 

1.8xl0 7 

3 . 9xl0 9 


Planetary Gear 
Characteristic 

Sun 

Gear 

Planet 

Gears 

Ring 

Gear 

Carrier 

Number of Teeth 

27 

35 

99 

— 

Diametral Pitch* 

8.8571 

8.8571 

9.14286 

- 

Pressure Angle* 
(degrees) 

24.6 

24.6 

20. 19 

- 

Root Radius (in) 

1.39 

1.802 

5.594 

- 

Tooth Tip Radius (in) 

1.655 

2.067 

5.35 

- 

Face Width (in) 

1.375 

1.20 

1.00 

- 

Inertia (in-lb-s 2 ) 

0.026 

0.016 

0.0 

0.20 


Note* The variation in diametral pitch and pressure angle cannot be 

accurately modeled, as the program assumes equal pressure angles 
at the sun-planet and ring-planet meshes. 
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Table 3: Two Stage Test Cas e Parameter Variations. 


CASE NUMBER 

ATTACHED INERTIAS 
(in-lbs-s* ) 

SHAFT STIFFNESSES 
(ln-lb/rad.) 


INPUT 

OUTPUT 

INPUT 

OUTPUT 

CONNECTING 

1 

10 

15 

10 

10 

10 

2 

100 

10 

10 

10 

10 

3 

15 

10 

10 

10 

10 

4 

10 

10 

10 

10 

lxlO 9 

5 

10 

10 

0 

10 

10 

6 

10 

10 

10 

10 

lxlO 5 

7 

10 

10 

lxl 0 6 

10 

10 
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Table 4: Results of Tvo Stage Test Case Parametric Variations 


DISPLACEMENTS, IN. 


CASE 

STAGE 

INPUT 

OUTPUT 

SUN 

CARRIER 

PLANET 

NUMBER 


SHAFT 

SHAFT 

GEAR 

5 

4 



Xl<J 

xlO 2 


xlO° 

xicr 

1 

1 

.8017 


.6941 

.7872 

.2067 


2 

— 

.8345 

.4274 

.503 

. 1292 

rs 

& 

1 

.08018 

— 

.06955 

.07985 

.02071 


2 

— 

.8018 

.6408 

.7528 

. 1937 

3 

1 

. 5345 

- 

.463 

.526 

. 1379 


2 

— 

.6018 

.6408 

.7534 

. 1937 

4 

1 

.8003 

— 

.6899 

.7629 

.2051 


2 

1 — 

.8003 

.639 

.7629 

.1928 

5 

1 

.8004 

— 

.616E-4 

. 1095E-2 

. 1974E-4 


2 

— 

.8004 

.638 

.7466 

.1925 

6 

1 

.8017 

— 

.6935 

.777 

.2065 


2 

— 

.8017 

.6414 

.764 

. 1938 

7 

1 

19.889 

— 

2755.6 

9612 

1247.9 


2 

- 

13.364 

59.4 

231.8 

26.1 


NOTE: These displacements are at the end of the first boundary 

condition iteration. 
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Stage 1 Stage 2 


Figure 1 : Two Stage System Diagram 


out 

a 

^outsh 
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Figure 2: Dynamic Model for Epicyclic Gear Stage 




(continued) 


Figure 3 : General Program Flowchart 
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(continued) 



(continued) 


Figure 3: General Program Flowchart (continued) 
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(continued) 
* 


Loop 


Figure 


for 10 time cycles 


1 


Loop for time step 


>i 



□g j 


3: General Program Flowchart (continued) 
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TCI 


"| Stage 1 


TC2 


TC2 


Stage 2 


0.0 


0.025 
time, sec. 


0.05 


TC1/TC2 = 0.05/0.025 = 2.0 


Number of Boundary Condition Cycles 
for Steady State - 2 

TCI = Time Cycle for stage 1 
TC2 = Time Cycle for stage 2 


-TC2 


-*+*■ 


■ TC2- 


Stage 2 


\ 


TC1-4--TC1 TCl-^TCl^*- 


1 Stage 1 


i t 


I ♦ 1 — ♦ 1 — H i- 

0.0 0. 01124 0.0264 

time, sec. 

TC2/TC1 = 0.0264/ .00562 - 4.7 


i 

-t — I — H — 


10 x TC2 - 0.264 sec. 
47 x TCI - 0.264 sec. 


Number of Boundary Condition 
Cycles for Steady State - 47 


Figure 4: Two Stage Time Cycle /Boundary Condition Diagram 
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Sun Gear Velocity, 



Figure 5 : Sun Gear Displacement and Velocity Boundary 

Conditions vs. Number of Iterations 
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Figure 8: Velocity Boundary Conditions vs. 

Number of Iterations 



SINGLE 

STAGE 

CHECK CASE; OH-58 




I 
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1 . 
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24.6 



3 

6. 

27. 

35. 
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1.0 

5 
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.026 

.2000 

.00 

.0160 

5 
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2 

150. 

.01 
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1 

651. 

1 . 





0 - 1 . 


Figure 9a: Single Stage OH-58 Test Case Input 
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Figure 10: Two Stage OH-58 Input Data Set 
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Figure lib: Dynamic Loads for OH-58 First Stage 

Ring-Planet Mesh 
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Figure 13: Procedure to Locate Source of Instability 
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Stage 1 


Stage 2 


Figure 14: Simple Two Stage System Schematic 
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